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Abstract. Second-order divergence- form operators with stationary random 
coefficients homogenize over large scales. We investigate the effect of certain 
perturbations of the medium on the homogenized coefficients. The perturbations 
that we consider are rare at the local level, but when occurring, have an effect 
of the same order of magnitude as the initial medium itself. The main result 
of the paper is a first-order expansion of the homogenized coefficients, as a 
function of the perturbation parameter. 
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1. Introduction 

Consider a second-order divergence-form operator whose coefficients are random. 
The randomness of the coefficients model the small-scale irregularity of a medium. 
If the distribution of the coefficients is stationary and ergodic (and if some ellipticity 
condition holds), then this random operator can, over large scales, be replaced by 
an effective operator with constant homogenized coefficients. 

The aim of this paper is to study the effect of certain perturbations of the medium 
on the homogenized coefficients. The perturbations that we consider are small only 
in the sense that locally, the medium is perturbed with small probability; but where 
a perturbation occurs, the change is of the same order of magnitude as the medium 
itself. This type of perturbation may be called a Bernoulli perturbation. Our main 
purpose is to prove a first-order expansion of the homogenized coefficients, as a 
function of the perturbation parameter, under conditions of short-range correlations 
and uniform ellipticity of the medium. 

There are at least two important motivations behind this problem. The first 
concerns the numerical approximation of the homogenized coefficients. Although 
several techniques for doing so have been identified and analysed [EGMN12], it 
remains a computationally expensive task, even in low dimensions. It has thus 
been proposed in [ALlOa] to study more efficient techniques that would apply to 
weakly random media, that is, random perturbations of a periodic environment. 
First-order expansions (as a function of the perturbation parameter) have been 
proved in [ALlOb], but only for specific types of perturbations that do not include 
Bernoulli perturbations. Yet, Bernoulli perturbations are arguably the most natural 
modelling assumption for typical disordered media like composite materials (as 
an example, see [ALII] for a cross-section of a composite material used in the 
aeronautics industry). In [ALII, ALlOb], conjectures are formulated concerning the 
expansion of the homogenized coefficients for such perturbations, which are backed 
by a formal derivation and numerical evidence. 

A second motivation is related to percolation. For this model, each edge of Z d 
is independently removed with probability p, and kept otherwise. There exists a 
critical probability p c € (0, 1) such that the remaining graph has a unique infinite 
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connected component if p < p c , and has only finite connected components if p > p c . 
Over the last years, the understanding of two-dimensional percolation close to 
criticality progressed tremendously, in particular through the rigorous derivation of 
the values of several critical exponents (see for instance [SW01]). Yet, to the best of 
my knowledge, it is not known (for any d ^ 2) how the homogenized conductivity 
of the percolation cluster behaves as p approaches p c from below (see [Hu09] for a 
nice review of the problem). Understanding the effect of small perturbations of the 
medium on the homogenized conductivity seems to be a first necessary step towards 
the resolution of this problem. 

Of the two motivations above, the first is formulated for differential operators, 
while the second is inherently discrete. In what follows, we focus on a discrete model. 
Apart from this difference, the results presented here give a proof of the first-order 
expansion of the homogenized coefficients conjectured in [ALlOa, ALII, ALlOb]. It 
extends it in the sense that the non-perturbed environment considered here need 
not be periodic, and that the expansion is obtained around every value of the 
perturbation parameter. We work under a condition of uniform ellipticity, which is 
obviously not satisfied in the case of percolation. The present paper will hopefully 
lay the basis for an extension to the case of percolation, as well as to higher-order 
expansions. 

Let us now give a brief survey of related works, outside of the previously mentioned 
[ALlOa, ALII, ALlOb]. 

In some cases, asymptotic formulas for the homogenized coefficients in the regime 
of high dilution of the perturbations have been known and used by physicists 
since the 19th century, and are variously known as the Clausius-Mossotti, Lorentz- 
Lorenz, Maxwell, or Rayleigh formulas. The setting is that of a continuous and 
homogeneous medium with highly diluted spherical inclusions. In this context, 
the general expression for the first-order correction that will be obtained here in 
terms of correctors becomes an explicit formula [Ma, Pa95]. Vanishingly small 
spherical inclusions arranged periodically along a fixed square or cubic lattice within 
a homogeneous medium have been investigated with precision in [St], where a high- 
order expansion is obtained. Surprisingly, in dimension 3, the expansion involves 
fractional powers of the volume fraction p of the small spherical inclusions, the 
smallest fractional power with a non-zero coefficient being p 4+1 / 3 (see [St, (60) and 
(64)]). In particular, the homogenized coefficients can be differentiated four times 
with respect to p, but no more. Note that in [St], every periodic cell is perturbed 
by the insertion of a small inclusion, while our present approach would cover the 
different situation where each periodic cell is perturbed by some fixed inclusion, with 
small probability. Yet, it is to be expected that in our context also, the homogenized 
coefficients are not infinitely differentiable functions of the perturbation parameter. 

Clausius-Mossotti formulas have been the subject of an impressive amount of 
work in the physics literature, of which we simply quote [ZB77], where results of [St] 
are extended to higher-order corrections. In the mathematics literature, the most 
important result to date on this problem (that I know of) is certainly [A113]. Under 
some conditions, the Clausius-Mossotti formula is proved (although the results 
obtained there do not directly relate to the homogenized coefficients), with an error 
bound of order p 3 ^ 2 , where p is the volume fraction of the perturbation. As we 
will see below, in our context and for some cases, we can show that the first-order 
expansion of the homogenized coefficients holds with an error term of order o(p 2 ~ ri ), 
for every rj > 0. Earlier investigations include [Ko89, BM01]. 

The problematic considered in [BLL07] is similar in spirit to the present one. In 
this work, a different type of perturbation, based on random deformations of the 
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geometry of the medium, has been investigated, and a first-order expansion has 
been obtained, see [BLL07, Theorem 3.2]. 

The problem of showing the regularity of diffusion coefficients has also been 
explored for others models, and in particular for that of a tagged particle in the 
simple symmetric exclusion process. It is proved in [LOVOI] that the effective 
diffusivity for this model is an infinitely differentiable function of the density of 
particles (see also [Be02, LOV04, Bc05, Su05, Na05, Na06, Na07] for generalizations). 
The approach followed there relies on a particular duality structure of the process. 
In our present context, this duality has been investigated in [KO05, CK08] under 
the additional assumption that the random coefficients take only two possible values 
(although no regularity result on the diffusion coefficients was given there). In 
contrast, the approach of the present paper does not involve any duality structure, 
and (as a consequence) is not restricted to random coefficients taking only a finite 
number of possible values. 

In the next section, we introduce some definitions and notations, and then explain 
the approach taken up and the organization of the rest of the paper. 

2. Definitions and notations 

2.1. A reminder on homogenization. We view 1 d (d ^ 2) as a graph, with its 
usual nearest-neighbour structure, and write B for the set of (non-oriented) edges. 
We write x ~ y if x and y £ % d are neighbours. Let (w e ) eS B be i.i.d. random 
variables such that almost surely, 

C- ^ W e ^ C+, 

where < c_ < c + < +oo are constants, henceforth called the ellipticity constants. 
We call ui — (w e ) e6 B an environment, and write fi = [c_,c+] B for the set of 
environments. We call w e the conductance of the edge e. We write P for the law of ui, 
and E for the associated expectation. Note that 7L d acts on f2 by the translations 
(0*)*ez<* defined by (9 Z u) XiV 

Let (ei, . . . , e<j) be the canonical basis of M. d . For every x £ Z d , let A(x,oj) = 
A u (x) be the d-dimensional diagonal matrix diag(u>a; x _|_ ei , . . . , uj x x + ed )- We may 
keep the variable ui implicit in the notation. The operator whose homogenization 
properties are investigated here is —V* ■ AV, which acts on functions / : Z d — > R by 

-V* • AVf(x) = "*,y(f(y) ~ /(*)) (* € 

Here and below, we write V/ to denote the forward gradient, 

" f(x + ei )-f(x) 

V/(s) = : 

_ f( x + e d )-f(x) 

and for F — (Fx, . . . , Fd) : Z d — > R d , we write V* • F for the backward divergence, 

d 

V* • F(x) = ( F *( x ) ~ p i( x ~ e «)) • 

Homogenization refers to the fact that there exists a constant symmetric matrix 
^4hom such that the solution operator of —V* • A(-/e)V converges, as e — > 0, to 
the solution operator of the differential operator —V • AhomV. Moreover, the 
homogenized matrix ^4hom can be characterized in terms of a function called the 
corrector, which we now proceed to define. 

We say that a function ip : Z d x ft — > K is stationary if tp(x, ui) — ip(Q, 9 X lu). Let 
^ £ M. d be a vector which will be kept fixed throughout this paper. The corrector <f> 
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in the direction £ is the unique function: Z d X — > K such that V0 is stationary, 
E[V0] = 0, 0(0) = and 

(2.1) - V* • A{£ + V<P)(x,oj) = (a € Z d , P-a.e. w). 

There are several ways to define the homogenized matrix in terms of the corrector, 
but the most useful one for our present purpose is the property that 

(2.2) e.A ham e = E[f.A(e+v^)]. 

In the right-hand side above, we keep implicit the fact that the quantity under the 
expectation is evaluated at (0, oj). Proofs of existence and uniqueness of as stated 
above, and of the formula (2.2), can be found for instance in [Kii83, Theorem 3 and 
(3.17)], or in the continuous setting, in [PV81, Theorem 2]. 

2.2. Bernoulli perturbations. We now introduce Bernoulli perturbations. We 
give ourselves a second family )eeB °f i-i-d- random variables such that almost 
surely, 

c_ lu^ c+. 

For every p G [0, 1], we also give ourselves a family (b e )eeB of independent Bernoulli 
random variables of parameter p, independent of (uj e ,uie)eEV,- We now understand 
that P denotes the joint law of lj := (ui e , u?', (&e ) P £[o,i])eeB) so that for instance, 
we have p = ¥[bi pS> = 1]. We write f2 for the set of possible values taken by lj. The 
group Z d acts on by translations, and we keep writing (0x) x ez d to denote this 
action. 

We let to e p ' = (1 — be) u) e + bi P ^ u>i . Note that the notation is consistent when 
p = 1, and that if we denote the law of ui^ by v^ p \ then = (1 — p)v^ + pv^ . 
With each p £ [0, 1] is thus associated the perturbed environment uj^ — (o4 P ^) e £B- 
The environment a/ p ) shares the same properties as the environment to, and thus we 
can define A^\A ( ^ m and <f>&) in the same way as A, Ah om , and 4> respectively, but 
with us replaced by ui( p \ Throughout this article, for notational convenience, we will 
replace the exponent (°) by simply °, so that for instance, A° — A^°\ A^ om = , 
cf)° = (f>(°\ and so on (we think of A° , 0°, etc. as functions of cj, which makes them 
formally different from A, 4>, etc. introduced before). 

Our main goal is to show that the homogenized matrix Ai?J , as a function of p, 
is differentiable, and to find an explicit formula for the derivative. Heuristically, 
one may expect that a linear approximation in (2.2) gives the correct first-order 
approximation, that is, as p tends to 0, 

(2-3) £ • 41C = £ ' ^U + pEeK • A e (£ + V<f ) - £ • A°(£ + V0°)] + °(p), 

where A e and <j) e are for the environment perturbed at the edge e what A and 4> are 
for the unperturbed environment. In the expectation above and throughout this 
paper, it is kept implicit that the functions are evaluated at (0, of). Our aim is to 
justify this heuristic. Note that the meaning of the sum in (2.3) is not clear, since it 
is not absolutely summable. 

Organization of the paper. In section 3, we recall classical results on the decay 
at infinity of the Green function and its gradients. A difficulty with the formula 
(2.2) is that the corrector is a very non-local function. It turns out convenient to 
introduce a localized version of the corrector, obtained by adding a zero-order term 
in equation (2.1). In section 4, we give a simple criterion for the existence and 
uniqueness of solutions of elliptic equations with a zero-order term. In section 5, 
we recall results quantifying the accuracy of the approximation by this localized 
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corrector. Although the (localized) corrector depends non-linearly on the values 
taken by the environment u)( p ', we will show in section 6 that it is, at first order, 
close to its linear approximation in the limit of small p. After having derived several 
useful convergence results related to stationary elliptic equations in section 7, we 
state and prove Theorem 8.1 in section 8, which is a rigorous version of (2.3). We 
show in section 9 that localization by a zero-order term in (2.1) can be replaced 
by periodization of the medium, thus providing us with alternative descriptions of 
the derivative of A-j*£ m at p = 0. The goal of the last two sections is to generalize 
Theorem 8.1 in two directions. In section 10, we give the asymptotic expansion of 
^hom as P tends to 0, for every p e [0, 1]. At the price of an additional assumption, 
we prove in section 11 that the error term in the expansion is o{p 2 ~ v ) for every 
7] > 0. 

Although we work throughout under the assumption that the medium is made 
of i.i.d. random variables, the method can be generalized to some (stationary) 
correlated environments, for instance those with a finite range of dependence. For 
this reason, we will not use any property specific to the i.i.d. case (as e.g. the fact 
that the homogenized matrices are multiples of the identity). 

Notation. We define aVt = max(a, b). We write | • | for the L 2 norm on M. d . This 
norm induces an operator norm ondxrf matrices, which we also denote by | • |. For 
a function of several variables, e.g. / : K x % d x Z d — > E, we write V2/ and V3/ to 
denote the (forward) gradient with respect to the second and to the third variable 
respectively. For i,j G {2,3}, we write ViVjf to denote the d x d matrix whose 
columns are Vj.Fi, . . . , V l F d , where [Fx,..., F d ] T = V f. Note that if v G R d is a 
fixed vector, and h = Vjf ■ v, then Vj/i = (VjV//)u. 

The value of a constant denoted by c or c may change from one occurrence to 
another, but depends only on the dimension and the ellipticity constants. 



3. Parabolic and elliptic regularity 

The aim of this section is to gather known results on the regularity of the heat 
kernel and Green function associated to divergence-form operators. We begin by 
defining the heat kernel. For every fixed uj £ f2 and x G 1 d , let io^it, x, y))(t,y)eR + xz d 
be the unique bounded function satisfying 

3 t ct{t, x, y) = -V* • A"(y)Vq"(t, x, y) ((i, y) G M+ x Z d ), 
F(0,x,y) = l y=x (yeZ d ), 

where we understand V* and V as acting on the y variable. Note that q u {t, x, y) = 
q^ (t,y,x). Let q*(t,x) = q^ (t,0, x), where lo* is the constant environment such 
that w* = 1 for every eel. 

Proposition 3.1 (pointwise control of the heat kernel [SZ97, De99]). There exist 
constants c, k > 0, a > such that for every u £ ft, t ^ 0, x G Z d and i G {2, 3}, 

(3.1) q u (t,0,x) ^ c q*(kt,x), 

q*{kt, x) 



(3.2) |V i g w (t > 0,s)|<c 



(1 Vt) a 



Proof. The first part appears in this form in [DD05, Section 4], and the proof given 
there relies on the closely related [De99, Proposition 3.4] (a similar statement was 
also obtained in [SZ97, Lemma 1.9]). It also appears that as soon as (3.1) holds for 
some value of k, it holds for every larger value as well (with a different c). This will 
also be clear from the proof of (3.2), to which we now turn. 
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Inequality (3.2) roughly corresponds to [SZ97, Theorem 1.31]. In order to provide 
the reader with a precise proof, we begin by recalling the following consequence of 
the parabolic Harnack inequality. For r > and x G Z d , let 

Q{x, r) = [0, 4r 2 ] x {z G 7L d : \z - x\ < 2r}, 
Q'(x, r) = [3r 2 , 4r 2 ] x {z G Z d : \z - x\ ^ r}. 
For a real function u defined on a set A, we write 

Osc it = sup u — inf u. 

A A A 

Proposition 3.2 (increase of oscillations [SZ97]). There exists A > 1 such that {or 
every oj G CI, r > and £ G Z d , j/u : R + xZ^R satisfies 

d t u = — V* • A(',w)Vti inQ(x,r), 

then 

Osc u ^ A Osc u. 

Q(x,r) Q'(x,r) 

The proof of Proposition 3.2 can be found in [SZ97, Lemma 1.30], or can be 
derived from the parabolic Harnack inequality given in [De99, Theorem 2.1] as in the 
proof of [DD05, Proposition 3.4] (see also [De99, p. 188] on the fact that numerical 
constants appearing in the definitions of Q and Q' are irrelevant, and can thus be 
chosen as convenient). 

We will also need the following result, which we borrow from [DD05, Section 4]. 
The expressions involved are less simple than what one might imagine a priori, due 
to lattice effects. 

Proposition 3.3 (explicit estimates on q*). Let 
T> t [x) — \x \ arsinh 




There exist constants c\ , C2 , fci, k 2 such that for every t ^ and x G 



■ exp(-T> klt (x)) ^ q*(t,x) ^ — — ^-T 7? exp(-2? fe2t (x)). 



(i v o d/2 v m " " ' "(IV ty/z 

Remark 3.4. The following alternative expression for T> t (x) is also given in [DD05]: 
V t {x) = max (s\x\ — t(cosh(s) — 1)) , 

which makes it transparent that 1 1-> T> t (x) is decreasing (in the wide sense). 

Let us see how to use these results to prove that (3.2) holds. We fix k such that 
(3.1) holds, and hi, k 2 as given by Proposition 3.3. We let fc 3 > be some parameter 
yet to be fixed, and begin by proving that (3.2) holds provided |x| < k 3 t. (This is 
the most relevant case, but lattice effects prevent us from giving a unified proof.) 

It follows from Proposition 3.3 that there exists constants c\,C2 such that if 
|o;| ^ 2kzkt, then 

V ; (1 Vt) d / 2 H V ' ; (1 Vt) d / 2 

Let x,t satisfy \x\ ^ k 3 t, and let L = [log 2 (i /2) / 2j . For L < (i.e. t < 2), there is 
nothing to prove since (3.1) holds. Otherwise, for every I G {0, . . . , L}, we consider 
the cylinders 

Qi 

Let y be a neighbour of x. Noting that (t,x) and (t,y) both belong to Qo and 
applying Proposition 3.2 iteratively, we obtain that 

|g w (t, 0, x) - q"(t, 0,y)\ < A~ L Osc<f (-,0, •) < \~ L supg"(-, 0, ■)• 

Ql 



= [t- 2 2l ,t] x {z e Z d : \z - x\ s? 2 1 }. 
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By the definition of L, the cylinder Ql is included in the cylinder 

Q = [t/2, t] x {z G 7L d : \z - x\ < Vt}. 

In view of (3.1) and (3.3), we thus get that 

\q u (t,0,x)-q u (t,Q,y)\^c\- L q*{k t,x), 

for every ko > large enough, uniformly over x and t satisfying |x| ^ k%t. This is 
indeed the inequality (3.2) by the definition of L and the fact that A > 1. 

There remains to justify (3.2) for |i| > k$t. Recall that we still have the freedom 
to choose &3 as convenient. By (3.1) and Proposition 3.3, we have 

(3.4) <f(t,0,x) <<f(fci,x) ^ v C /2 exp(-P fc2fct (x)). 

We now claim that choosing k 3 sufficiently large, we can ensure that for every 
k > 2kk 2 , 

(3.5) \x\ > k 3 t V k2kt {x) ^t + V u (x). 

Let us first see why (3.5) enables to conclude. It follows from (3.4) and (3.5) that 

q»(t,0,x) < v C f)d/2 exp (-P fa (s)) . 

By Proposition 3.3, this is smaller than 

ce q*(kt/ki, x), 

for every k > 2kk 2 and every a;, t satisfying |x| ^ fc 3 t. This clearly implies (3.2) for 
|x| ^ k 3 t. The proof is thus complete, provided we show (3.5). 
In order to do so, we first note that 

„ -i)-VP+HF-«-jf TO d "« w ' 

The right-hand side of (3.5) thus holds as soon as 

(3.6) \x\ arsinh — — ^ t + \x\ + \x\ arsinh . . , 

\k 2 kt J \zk2kt 

If k 3 is sufficiently large (and in particular ^ 1), then it is clear that |x| ^ k 3 t 
implies 

arsinh — - ^ 2 + arsinh ' ' 
which in turn implies (3.6). □ 

The bound (3.2) on the gradient of the heat kernel cannot be improved in general 
(see for instance [DD05, p. 364]). However, upper bounds matching the homogeneous 
case can be recovered if one takes averages over the randomness of the stationary 
field of conductances. For simplicity, we write q^ p \t, x, y) instead of g w<P) (t, x, y). 
The following result is due to [CNOOa, Theorem 1.4] and [DD05, Theorem 1.1]. 

Proposition 3.5 (averaged control of the heat kernel [CNOOa, DD05]). There exist 
constants c,k > such that for every p £ [0, 1], t > x G Z d and i,j G {2, 3}, 




(3.7) (E[\V iq ^(t,0,x)\ 2 }) 



1/2 q*(kt,x) 
< c — , , 

Vi vt 



(3.8) E[|V 4 V,V p) (t,0,x)|] S =: c ffl^ ) 
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For every u> £ f2, p, > and x, y € Z d , we define the Green function 

f+oo 

G"Jx,y) = / e-^q u (t,x,y) dt. 



Jo 

Note that the function (G^(x,y)) yeIi d satisfies 

»G»(x, y) - V* ■ A(y, ^VG^x, y) = l x = y (y e Z d ), 

where V* and V are understood as acting on the y variable. Note also that the 
Green function is symmetric, that is, G^(x,y) — G^(y,x). The estimates obtained 
on the heat kernel transfer into estimates on the Green function by integration. 

Proposition 3.6 (pointwise control of the Green function). There exist c,c > 
0, a > such that for every /i G (0, 1/2], oj <E Q, x e Z d and i e {1, 2}, 



(3.9) \G"(p,x)\^c 



log^- 1 ) e-^l ifd=2, 

1 e -S>/P|x| j/rf^ 3j 



(iv N) 



rf-2 



(3.10) |V 4 G^(0,x)| *S (lv| ^ )d _ 2+a e- E vW. 

Proof. We begin by justifying (3.9) for d ^ 3. In view of (3.1), it suffices to bound 

r+oo <.\x\ r+oo 

(3.11) / e-^p*(t,x) dt = / e-^p*(t,x) dt+ e^^p*(t,x) dt. 

JO JO J\x\ 

Using (3.3), we can bound the second integral in the right-hand side above, up to a 
constant, by 

f+OO p —ixt 

(3.12) / w/9 e~ N /C2t dt, 



(lVt)^- 

which, by a change of variables, can be bounded by 



r+oo -/j,\xfu 



Moreover, 



/■(v^M) „-l/2c 2 « r+oo -l/c 2 « 

^ e -VPN/2c 2 / ^^ du + e -VP\*\ —rp^&u 

Jo u d / 2 J(^l)- 1 u 

^ (7( e -v / ^l ;r l/ 2c 2 + e-x/^l^l). 

We have thus obtained that for d ^ 3, there exists c, c > such that 

'|*| Fl 

The first integral in the right-hand side of (3.11) can be bounded more easily. Indeed, 
it follows from Proposition 3.3 that for t < \x\, we have p*(t,x) ^ e~ c \ x \ for some 
c > 0, and thus 

"M 

e-"V(t,x) dt < \x\e- £lxl . 

o 

To summarize, we have thus shown that 

e-"V(*,x) dt < ri5=2 e " eVP|a:| - 
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In order to conclude, it then suffices to see that the Green function is bounded 
uniformly over fi. But this is clear since Jp*{t,0)dt is finite. 

When d — 2, one needs to be more careful when considering the integral appearing 
in (3.12). It is straightforward to check that this integral is bounded by a constant 
times log(^ _1 ), uniformly over x. If y^l^l ^ 1, then there is nothing more to prove. 
Else, we split the integral along 



-fit 



M/vW 1 Vt 



-\x\ 2 /c 2 t ^ 



The first integral is bounded by e ^ / C2 , and can thus be discarded. For the second 
one, we obtain the bound 



t 



dt sC 21og(^~ 1 )e~ 



/MM/C2 



while the third integral can be bounded by 

f+OO 



/MM/2 



e -Att/2 e -^Tfi\x\ 

dt ^ 



'MR 



and the claim is thus proved. 

Inequality (3.10) is obtained from (3.2), following the same reasoning as for the 
case d ^ 3 of the proof of (3.9). □ 



We write for G" (p) . 

Proposition 3.7 (averaged control of the Green function). There exists c, c > 
such that for every p £ [0, 1), p > 0, x £ Z d and i,j £ {1, 2}, 



(3.13) 



(E[|ViGW(0,x)| s 



1/2 



(iv M) 



d-l 



(3.14) 



E[| 



V«V^)(0,x)|]<^^_e- 



Proof. Recall that Minkowski's integral inequality ensures that for any positive 
measurable function /, 



E 



f(s) ds 



Applying this observation to 

E[\ViG<*\0,x)\ 2 ] sSE 



E[/( s ) 2 ] 1/2 ds 



-^\V iq ^(t,0,x)\ dt 



we get that 



1/2 



+ 00 

: / e - "' E 
o 



\V iq M(t,0,xtf 



1/2 



dt. 



The rest of the proof of (3.13) is the same as that of (3.9), proceeding by integration 
of the estimate obtained in (3.7). The proof of (3.14) is identical, integrating 
estimate (3.8). □ 
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4. Solving regularized elliptic equations 

We say that a function / : Z d — > M is of sub-exponential growth if for every 5 > 0, 
we have |/(:r)| = 0(e*\ x \) as I a; I tends to infinity. 

Theorem 4.1 (Existence and uniqueness of solutions). Let fi > 0, and f : Z d — » K 
&e 0/ sub- exponential growth. For every lu G SI, £/iere exists a unique function 
X : Z" 1 -) I 0/ sub- exponential growth satisfying 

(4.1) «C - V ■ A u Vx = / (mZ d ). 

Moreover, x is given by 

(4-2) X (x)= £G-(x,y) /(y) (x G Z d ). 

Proof. Since (3.9) ensures exponential decay of the Green function, it is easy to see 
that the function \ written in (4.2) is well-defined, of sub-exponential growth, and 
satisfies equation (4.1). There remains to show uniqueness, and for this, it suffices 
to consider the case / — 0. Let B n = {— n, . . . , n} d be the box of size n. We have 

(4.3) Y, ^ M x ) ~ V * ' A u {x)V X {x)) = 0. 
xeB n 

We want to perform an integration by parts. The difference between 

- ]T x(z)V* • A"(x)V X (z) 

x£B n 

and 

(4.4) J2 Vx(x)-A u (x)V X (x) 



xeB„ 



is bounded by 



J2 \x(x)\ \A"(y)V X (y)\ + \x(y)\ \A u (x)V X (x)\, 



x<£B n ,yeB n+1 
x~y 

which, up to a constant, is bounded by 

xeB n+2 \B n - 2 

Hence, combining this estimate with (4.3) and the fact that the term in (4.4) is 
positive, we obtain 

x£B n ieB„ +2 \B n _ 2 

for some constant c > 0. In a more condensed form, if we write 

we have thus obtained that /1 u„ ^ c(u n+ 2 — u n -2)- Since u n is increasing, we have 
shown that 

U n +2 > (l + "J Un-2- 

If w„ is not identically zero, then it must grow exponentially fast. But this would 
contradict the assumption that x is °f sub-exponential growth. Hence u n is identically 
zero, and so is x- D 
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5. Approximating the homogenized matrix 
Let fx > 0. For every p G [0, 1] and w G O, we let ^ p) (-, w) : Z d 



be the 



unique function of sub-exponential growth satisfying 
(5.1) ^)-V*-AW(e + V^)=0. 

We write 6° for 6}?' ■ The introduction of the function 



6° for . The introduction of the function </>[f is interesting for our 

purpose as an approximation of <f>^ p \ The former is more localized than in the 

(v) 

sense that, roughly speaking, the value of </>ji at a point depends only on the values 
of the conductances in a box of size /i -1 / 2 . 

A standard energy estimate shows that E[|V0^f | 2 ] is bounded uniformly over p, 
and p (we recall that whenever we write an expression like E[| V<f)p | 2 ], we understand 
that the function under the expectation is evaluated at (0,w)). We will need the 
following much more precise information on the corrector, which is borrowed from 
[GOll, Proposition 2.1]. 

Theorem 5.1 (finite moments of the corrector, [GOll]). For every k G N ; there 
exist constants c and r ^ such that for every \x G (0, 1/2] and p G [0, 1] , 

log^- 1 ) ifd = 2, 
1 ifd^Z. 



E 



1^1* 



^ c 



It is crucial for our subsequent arguments to have a good quantitative estimate 
on the difference between the homogenized matrix and the approximation obtained 



through 



The following result is what we need, and follows from [G012, 



Theorem 1] (or also from [GOll] and [Moll, Proposition 9.1 with k = 1]). 

Theorem 5.2 (approximation of Ahom based on [GOll, Moll, G012]). There 
exist constants c and r ^ such that for every fi G (0, 1/2] and p G [0, 1] , 



A ip) £ 



•E[£.A<p>(£ + V0W)] 



«S c 



plog^p- 1 ) ifd = 2, 
p if d ^ 3. 



6. Linear approximation of the corrector 

Each edge e can be written uniquely as (z, z + &i) for some z G Z d and 1 sC i d. 
We write e = z. We define C e : Z J x SI -) (the set of d X d matrices) by 
C e (x, of) = (the matrix whose elements are all 0) if x ^ e, and C e (e, lo) to be the 
diagonal matrix whose diagonal elements are all except the z th one, which is equal 
to wjp — uj e . As usual, we may keep the variable u implicit if clear from the context. 
We let A e = A° + C e . In words, the environment corresponding to A e is perturbed 
only at the edge e. More generally, for every E C B, we define 

(6.1) C E = J2° e and A E = A° + C E , 

e<EE 

G E for the Green function corresponding to the environment given by A E , and 
Eb) = {e G B : bi p) = 1}. Note that by definition, = A eM and = G e1p> . 
For every E C B, we let cj> E : Z d x O -> R be the unique function of sub-exponential 
growth such that 

(6.2) /<-V*.A B (£ + V^)=0, 
as given by Theorem 4.1, so that 

(6.3) 0» = J2 G >> V) V * ' AE (yK = E V ^>, y) • 
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where the integration by part in the last equality is justified since G decays expo- 
nentially fast, see (3.9). We have 4>^ — £<P> . We write <f>^ as a shorthand for (f>]f^ . 
We also define 

(6-4) #f = 0£ + E(#-^)' 

e£E 

(we will see shortly that this definition makes sense) , and write <^f for <^ <P) • The 
purpose of this section is to prove the following result, which roughly speaking, 
states that in the limit of diluted perturbations, the perturbed corrector is close to 
its linear approximation. 

Theorem 6.1 (linear approximation of the corrector). There exist constants c, 
ft > 0, 7 > such that for every \i 6 (0, 1/2] and p G [0, 1], 



E 



We start with several lemmas concerning the effect of perturbations on the 
corrector or the Green function. 

Lemma 6.2 (perturbation at one edge). Let 

(6-5) % = <t>%~ 

For every /i > and x £ Z d , we have 

%{x) = ~V 2 G^(x,e) -C e (e)(e + V^(e)). 

Proof. In view of equation (5.1) satisfied by <j>° and of the identity A e = A° + C e , 
we infer that 

- V* • A e (£ + V0°) = -V* • C e (C + V0°). 
This and equation (6.2) satisfied by <jf lead to 

(6.6) ii% - V* • = V* • C e (£ + V^°). 

The function on the right-hand side above is non-zero only at a finite number 
of points, so it is clearly of sub-exponential growth. The function ip^ is also of 
sub-exponential growth by (6.3) and (3.9). By Theorem 4.1, 

4>l{x) = X^foy) V*-C e (£ + V0°)(y) 
y 

y 

= -V 2 G^(x,e).C e (e)(£ + V0°(e)), 

where in the second equality, the integration by parts is justified since only a finite 
number of terms in the sum are non-zero. □ 

Remark 6.3. As was recalled in the beginning of section 5, the gradient of <f>° is 
square-integrable: E[|V</>°| 2 ] < oo. hence, by the ergodic theorem, 

(6-7) E tl V ^l] < 

It thus follows from Lemma 6.2 and the decay of the Green function (3.9) that for 
every x, fi^ix) decays exponentially fast with the distance from x to the edge e. In 
particular, for every E CM, the function <j>^ defined in (6.4) is always well-defined 
and of sub-exponential growth. 
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Lemma 6.4 (Perturbation at several edges). Let E CM, and let 

(6-8) $ = 

For every /i > and x £ Z d , we have 

4>ti*)=- E V 2 G^(x,e')-C™(e')V^(e'). 

Proof. If e G E, then using (6.6) and the fact that A E = A e + C E \ {e} , we have 

- V* ■ A E V% = V* ■ G e (£ + V0°) - V* • c E \^v%. 
Similarly, we have 

Ml - V* • A E (C + V0°) = -V* ■ G £ (£ + V«£°). 

Recalling that 

</V = 4% ~ 4>1 ~ E^V' 
that G B = J^esB anc ^ equation (6.2) satisfied by 4> E , we obtain 

— B 

By Remark 6.3, the function <j) is of sub-exponential growth. Moreover, arguing 
as in this remark, we see that the function in the right-hand side above is also of 
sub-exponential growth. By Theorem 4.1, we thus have 

eG-E y eZ d 

= -EE V 2 G E (x,y)-C E ^(y)V^(y) 

eEE yEZ d 

= E V 2 Gf(x,e')-C B \W ( e')V^(eO, 

where the integration by parts in the second equality is justified using the exponential 
decay of the Green function and the fact that y i— » \4>Ay)\ is of sub-exponential 
growth, see Lemma 6.2 and (6.7). □ 

Lemma 6.5 (averaged control of the perturbed Green function). For every £ € (0, 1], 
there exists c > such that for every jj, 6 (0, 1/2] and e, ef € B, 



E 



|V 1 V 2 G^(e',e)| 1 /C 



c 



(IV |e-e'|) d ' 



Proof. As a first step, we show that the lemma is true if the perturbed Green 
function G^ is replaced by the non-perturbed one G° . 
We know from (3.14) that 

Moreover, by (3.10), the quantity V2G^(x,y) is bounded uniformly over u e fl, 
x, y G 1 d and [i > 0. In particular, the random variable 

|ViV 2 G°(e',e)| 

is bounded. Hence, the lemma is true if G^ is replaced by G°. 
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Let us write = — G°. In order to prove the lemma, it now suffices to 
argue that 



E 



|ViV 2 G„(e',e)| 



(IV \e-e'\) d ' 



Let x E Z d . The Green function G°{x, •) satisfies 

vG° fl (x,-)-V*-A°VG° fi (x,-) = l x . 

The same equation holds for G e (x, •) provided A° is replaced by A e = A° 
hence 

fiG^x, •) - V* • A°VGl(x, ■) = lx + V* • C'VGfa, ■). 
Combining the two previous identities, we obtain that 

M G>, •) - V* • A°VG>, ■) = V* • rVG^, •)• 
By Theorem 4.1, we are led to 

K&v) = E G ^> Z ) v*'C e vG;( v )( z ) 

zSZ d 



£ V 2 G°(y,z) ■C e VG;( a; ,-)W 



= V 2 G°( 2 /,e) -G e (e)V 2 G^(x,e), 
and by symmetry of the Green functions, we also have 

(6.9) G>,y)=V 2 G°(x,e) • G e (e)V 2 G^, e). 
It thus follows that 

(6.10) ViV 2 G;(a ; ,y) = G e (e)(V 1 V 2 G^(2/,e)) (ViV 2 G^(»,e)) 
which we need to evaluate for x — e' and y — e. Since, by (3.10), 

|G e (e)VxV 2 G£(e,e)| 
is bounded uniformly over fj, and w, it suffices to see that 



E 



|ViV 2 G°(e',e)| 



(1 V |e-e'|) d ' 

but this was obtained during the first step of this proof, so we are done. 
Proof of Theorem 6.1. Note first that 



E 



where as usual, we write 



E 



(p) 



v^ p) (o) 



. By Lemma 6.4, we have 

(p)\ 



E (v 1 V 2 G(f)(0,e')) C EW \W(f')Vr,(e'), 

e^e'e,E(p) 

so we have the bound 

|V^ p) (0)| < c + J] |ViV 2 GW(0,e')| |V^(e')|. 

By Lemma 6.2, the right-hand side above is bounded, up to a constant, by 
(6.11) ]T iViVaG^J^e'JI |V x V 2 G£(e', e)| |£ + V0°(e)|. 



Let 7 > 1/2, and write E^ 



C e 



□ 



..,„ - fiW n S^-, , where we recall that £L = 
{—ft, . . . ,n} d . We will argue that in the above sum, the only terms which truly 
contribute to the sum are those for which e, e' € Ej?\ In order to do so, it is 
convenient to introduce 7 such that 1/2 < 7 < 7, and to define Eff' — E^ n B^-=,. 
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Recall that E[|V</>° (e)|] E[|V0° (e)| 2 ] 1/2 does not depend on e and is bounded 
uniformly over \i > 0. We also have from (3.10) that for every E CM, 



(6.12) 



\V 2 G^(x,y)\^ce- 5 ^\y- x \, 



with c > 0. Hence, up to a constant, the expectation 
(6.13) 



E 



E E |ViV 2 G|f)(0,e')| IVxVsG^e'.e)! |£ + W°(e)| 



is bounded by 

Summation over z gives a contribution of order ^T d / 2 , while the remaining summa- 
tion over z' decays as "' for every 7' < 7 — 1/2. We have thus justified that 

the difference between the expectation of (6.11) and (6.13) is 0(e~ M 1 ) for every 
7' < 7 — 1/2. We can thus focus on studying 



E 



E 

'e-E< p) ,ee_E<P),e^e 



|V 1 V 2 GW(0,e')| |ViV a G^(e',e3| |£ + V0°(e)| 



A similar reasoning enables to show that in the above sum, the terms for which 
can be discarded. Indeed, the contribution of those terms can be bounded, 
up to a constant, by 

(6.14) E 



E 

z'eB ^ z& d \B 



-c^fp,\z'\ -c^fp\z~z'\ 



Since 7 < 7, the distance between z and z' in any term of this double sum is always 
greater than /i~ 7 /2. Hence, we obtain an upper bound for the inner sum in (6.14) if 

we sum over all z such that \z — z'\ ^ /i~ 7 /2. This gives a sum of order 0(e~ M 7 ) 
for some 7' > 0, and hence a similar bound holds for the total sum in (6.14). We 
have thus argued that it suffices to study 



E 



E |VxV 2 GW(0,e')| |ViV a G£(e',e)| |f + V0°(e)| 



Finally (in order to keep light notations), arguing as in the first step, we see that 
we can as well consider the sum as ranging over all e,e' e E^ , e^e'. 
We now decompose the expectation into 



(6.15) E 



E |V 1 V 2 GW(0,e')| \^ 1 V 2 G;(e',e)\ |£ + V^(e)| 



= E E £(e,e',E)F{EW=E], 

E e=£e'£E 

where in the first sum, E ranges over all possible subsets of B„—y, and where we 
write 



f(e,e',£;)=E[|V 1 V 2 G^(0,e')| |ViV 2 G«(e', e)| |£ + V0°(e) 



fiJW = E 
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For any £ <E (0, 1), we have by Holder's inequality 



(6.16) £(e,e',E)^E 



(|V 1 V 2 GW(0,e')| |V 1 V 2 G^(e',e)| 

\t + v<f>°M 1/{1 - c) E M=E 



r 






A 1 



E 



i-C 



Since |£ + V0°(e)| does not depend on E^ p \ the last expectation is simply 

i-C 



E 



\t+v<t>°M 1/{1 ~° 



which is bounded by a power of log(^ _1 ) by Theorem 5.1. 

For the first expectation in the right-hand side of (6.16), we use the pointwise 
bound (3.10) on the gradient of the Green function to get 



E 



(|ViV a GW(0,e')| |V 1 V 2 G^(e',e)f 



r 


= E 




h 1 



^ (IV |e'|) rf - 2 + tt E 

r 

E 



V 1 V 2 G^(e',e)| 1 ^ = E 

V 1 V 2 G^(c',e)| 1 /C 



(1 V \e'\) d - 2 + a 
and Lemma 6.5 ensures that the latter is bounded by 



(IV |e'|) d - 2 +« (IV |e — e'|)C d " 

It thus follows that the right-hand side of (6.15) is bounded by some power of 
log(^t _1 ) times 



E 



^ (1 V \e'\) d - 2 + a (1 



(X v |e'|)d-2+a (i v \e-e'\)< d 



* 2 S n 



(1 V \d\) d - 2 + a (1 V le-e'l)^' 



To sum up, we have shown that there exists 7' > such that for every £ < 1, 



E 



< cp 2 log- (v- 1 ) Yl a 



1 



(1 V |e'|) d - 2 +« (1 V |e-e'|)f d 



for some exponent r ^ 0. Up to a multiplicative constant, the sum over e, e' above 
can be bounded by the integral 



/ 



dx dy 



|d-2+ Q | 2/ _ a ,|Cd 



'|x|</i-T,|j— a:K 2 /i — 1 \ x \ 



dx dy 



-2+a\ y _ x \(d 



Until now, the parameters 7 > 1/2 and £ < 1 were arbitrary. We can choose them 
in such a way that 

j(d - (d + 2 - a) < 1, 
and this finishes the proof of Theorem 6.1. □ 
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7. Solving stationary elliptic equations 

Recall that we say that a function if> : Z d x £2 —> R is stationary if ip(x,uf) = 
ijj(0,8 x w), where (0 x ) xe %d denotes the action of Z d on £2. The stationary function 
ip is thus fully characterized by the knowledge of the function u — > "0(0,0;), with 
which it can be identified. 

The first purpose of this section is to recall existence and uniqueness results for 
elliptic equations. Contrary to what was done in section 4, the equations considered 
here do not contain a zero-order regularizing parameter. On the other hand, we 
will work with the extra assumption of stationarity. Because of the identification 
mentioned above, we will focus on functions defined on £2 (instead of Z d x £2). We 
may keep this identification implicit, writing for instance A° instead of A°(0, •). 

For / : £2 — ► R, we define the forward gradient Df : £2 —> M. d as 

" f(o ei w) - f(u>) 

Df{io) = : 
and for F — {F\, . . . , Fd) : £2 — > R d , we write £>* ■ _F for the backward divergence, 

d 

D* ■ F{u) = J2 (Fi(u) - Fi(e_ ei w)) . 

i=l 

We write ||F|| 2 for the L 2 norm of F, that is, \\F\\% = E [ F il We let L%( M be 

the space of functions F : £2 —> R d such that ||F||2 is finite, and Ly be the closure 
in L 2 (£2) of {£)/, / : £2 R s.t. E[/ 2 ] < oo}. 

Theorem 7.1 (Existence and uniqueness of solutions). For every F £ i 2 (£2), there 
exists a unique \ € -^v such that 

(7.1) - D* ■ A° X = D*-F. 

Proof. We begin by showing the existence of solutions, following roughly the ar- 
guments of [KV86]. (Although the proof is by now standard, it is useful to recall 
its workings since we will need to extend it slightly later on.) The operator 
C := -D* ■ A°D is self-adjoint and positive on {/ : £2 -> R : E[/ 2 ] < oo}. The 
positivity comes from the observation that, for every square-integrable /, 

(7.2) E[/ £/] = E[Df ■ A°Df]. 

Hence, for every fj, > 0, there exists a square-integrable "f^ : £2 — > R such that 

(7.3) (// + = D* ■ F. 

Let f — D* ■ F. Since / is square-integrable, by the spectral theorem, there exists a 
measure e/ on R + such that for every bounded continuous function G : R+ — > R, 
one has 

(7.4) E[fG(C)(f)} = jG(X)de f (X). 

Note that for every square-integrable function g, one has 

E [/ g] = -E[F ■ Dg] < ||F|| 2 \\Dg\\ 2 , 
by the Cauchy-Schwarz inequality. Furthermore, by (7.2), 

\\Dg\\l < ^E[ ff £ 9 ], 



so that 



F\ 



E [/ 5 ] ^ ^ E[ff £ ff ] 



il/2 
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For every n £ N, let G„(A) = A 1 An. Taking g — G n {C){f) in the inequality above, 
and using (7.4), we get 

(A" 1 A n) de/(A) sc JS 2 - (^J A(A~ X A n) 2 de/(A)^ ' . 



By the monotone convergence theorem, this turns into 
(7.5) /I de/ (A)< l|F|1 ^ 



Equipped with this inequality, we can show that (£ )v E' A1 ) A1 >o is a Cauchy sequence in 
L 2 (f2). Indeed, note that by (7.2), we have 

(7.6) ||D* M - D9 v \\l < ^-E[(* M - ¥„) - *„)]. 

Since = (/i + £) -1 /, and using (7.4), we can rewrite the expectation in the 
right-hand side above as 

A^ + A^-^ + A)- 1 ) 2 dey(A) = j { ^ {V xn f +xy dey(A). 

Without loss of generality, we may assume that ji < v. In this case, the integrand 
above is bounded by 

Xu 2 _ 1 
X 2 ^ 2 - y 

which is integrable by (7.5). Now, for every fixed A, we have 

X{v-p) 2 Xv 2 



(^ + X) 2 (v + X) 2 ^ A 4 ' 

which tends to as /i, v — > 0. By the dominated convergence theorem, we thus obtain 
that the left-hand side of (7.6) tends to as /i, v — » 0, and thus that (-D5' M ) M >o is 
a Cauchy sequence in L 2 (Q). Let us write x f° r the limit. By definition, we have 

We now show that \ satisfies (7.1). One can check that, as a consequence of the 
estimate (7.5) (and using (7.4)), 

(7.7) ^ 0. 

For every square-integrable g : £2 — > K, we can write the weak formulation of (7.3): 

fjE\g + E[Dg ■ A°D^f ll ] = E[g D* ■ F}. 

Using (7.7) and the Cauchy-Schwarz inequality, we get that (jM[g tends to as 
/i tends to 0, and thus 

(7.8) E[Dg-A° X ]=E[gD*-F]. 

The left-hand side above is equal to —E[g D* ■ A°x]- Hence, the function D* ■ A°x + 
D* ■ F is orthogonal to every function in L 2 (fl). It is thus equal to zero, and that is 
to say that x satisfies (7.1). 

We now turn to uniqueness. By linearity, it suffices to show uniqueness for F = 0. 
Let x G Ly satisfy 

(7.9) -D*-A°x = 0- 

There exists a sequence of square-integrable /„ : —> M. such that Df n converges 
to x in L 2 (fl). By the weak formulation of (7.9), we have 

E[Df n ■ A° x ] = 0. 



EXPANSION OF HOMOGENIZED COEFFICIENTS 



19 



Passing to the limit, we get E[x • A°x] = 0, which implies that x — 0, and thus 
finishes the proof. □ 

We now proceed to derive several convergence results that will be useful for our 
subsequent reasoning. 

Proposition 7.2 (convergence of correctors). The function V0° (0, •) converges in 
L 2 (fl), as [i tends to 0, to the function V<^°(0, •), while for every eel, the function 
V0^(O, •) converges in L 2 (fl), as fi tends to 0, to a function that we write V</> e (0, •). 

Proof. The first part is classical (and is a minor adaptation of the proof of Theo- 
rem 7.1, choosing F — A°£). For the second part, in view of Lemma 6.2, it suffices 
to study the convergence of V 2 G^(x,y), for any fixed x,y £ Z d and u> € O. We 
define 



r+oo 

V 2 G u, {x 1 y) = / V 3 q u (t,x,y) dt. 
Jo 



Note that this definition makes sense even in dimension 2 by (3.2), although G u (x, y) 
itself does not. For the same reason, V2G"(x, y) is bounded uniformly over oj E O,. 
We get that 

r+oo 

V 2 G u (x,y)-V 2 G"(x,y) = / (1 - e^)V 3 ^(i, x, y) dt. 



This and (3.2) ensure that V 2 G^(x,y) converges to \7 2 G ul (x,y) as fi tends to 0, 
uniformly over w e 51. Using Lemma 6.2, we thus obtain the convergence of V<^ in 
L 2 (H), as desired. □ 

Remark 7.3. The proof shows that 

(7.10) V0 e (O,w) = V0°(O,o;) - (ViV 2 G e (0,e))C* e (e)(£ + V^°(e)). 
Proposition 7.4 (convergence of full correction). Let ip^ : ft — > M be defined by 

The gradient Dip^ converges in L 2 (fl), as \i tends to 0, to x € Ly i/ie unique 
solution of 

(7.11) - D* • A°x = D* ■ F, 
where F € £ 2 (H) is defined by 

(7.12) F = ^C e (£ + V<f )((),•). 

Remark 7.5. Note that ^ is well-defined, since 0® introduced in (6.4) is. Also, 
in the definition of F, note that the presence of the term C e ensures that every 
summand indexed by an edge e such that e ^ is actually equal to 0. 

Proof. Recall that, by definition, cj>^ satisfies 

^-V*-A e (£ + V^) = (inZ d ), 

and that A e = A° + C e . As a consequence, we get 

- V* • A°({ + V(j>*) = V* • C e (C + V^) (in Z d ). 

Using also the definition of <f>°, see (5.1), and recalling that we write <j> e = 0® — tj>°, 
we obtain that 

y4l - V* • A°VX = V* • C e (£ + V^) (in Z d ), 
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and thus 

(7.13) /i^-V*-A°VV^-^V*-C e (e + V^) (mZ d ). 
We let 

(7.14) iy^H^Tc^ + v^)^). 

e£B 

Note that is a stationary function, so we can identify it with the function 

-> R d 
lo i-> i^(0,w). 

With this identification in mind, we can rewrite (7.13) as 

(7.15) (n-D*-A°D)^ = D*-F M . 

By Proposition 7.2, the function Fu converges in L 2 (£2) to F defined in (7.12). We 
would like to argue as in the "existence" part of the proof of Theorem 7.1 to show 
that Di/jp converges to \ satisfying (7.11), but a difficulty arises since the right-hand 
side of (7.15) now depends on \i. Let ip^ € L 2 (fl) be the unique solution to 

(fi-D*- A°D) ^ = D* ■ F. 

By Theorem 7.1, we know that Dip^ converges in £ 2 (f2) to x € Ly the unique 
solution of (7.11). We now let ip„ = ip^ — ip^, so that 

(H-D*-A°D)'$ II = D*-(F M -F). 

In order to conclude, we need to show that Dip^ tends to in L 2 (f2). Let us write 
f li =D* ■ (Ffj, - F). In view of (7.2) and (7.4), we have 

PVv|| 2 < = ±-j de /M (A) < ±-j \ de u (X). 

By (7.5), the last integral is bounded by 

ll^-^lli 



and since F^ tends to F in L 2 (fl), this finishes the proof. □ 

8. First-order expansion around p = 

We are now ready to state and prove the first-order expansion around p = of 
the homogenized matrix Az\J . 

Theorem 8.1 (first-order expansion around p = 0). There exists a° S R such that, 
as p tends to 0, 

(8-1) C-4oU = ^-^ om C+pa?+ (p). 

Moreover, the coefficient can be defined as the limit of a\([i) as \x tends to 0, 
where aKn) is given by 

(8.2) a?(/x) = £ (Efc ■ A e (£ + V^)] - E[£ • A°(£ + V<£°)]) , 

eGB 

and where we recall that <p e ^ satisfies (6.2) wzi/i _E = {e}, that 0° satisfies (5.1) wrf/i 
p = 0, and i/iai i/ie functions under the expectations in (8.2) are understood to be 
evaluated at (0,a;). Alternative characterizations of a° are given in Propositions 9.1 
and 9.2 below. 
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Remark 8.2. The proof actually reveals that there exists an exponent fj > 
depending on the cllipticity constants such that (8.1) holds with o(p) replaced by 
o(p 1+7 '). Under an additional assumption, we show in Theorem 11.3 that the error 
term is actually o{p 2 ~ ri ) for every r\ > 0. 



Proof. For this proof, we fix 
(8.3) 



i" =P 



l+s 



for some e > 0. Note that in order to avoid heavy notations, this dependence 
between fi and p is kept implicit. For instance, and in view of Theorem 5.2, we may 
write without further comment that 



(8.4) 



hom> i> v*s 1 t pb 

For fj > given by Theorem 6.1, we have 
(8.5) 



(p->0). 



2-(l+ e )(l-,8) 



E[£ • A«(£ + V$f>)] - E[£- A<% + V#f>)] = O (p 
We fix the value of e > in such a way that 
(8.6) 2-(l + e)(l-0)>l. 

Using the definition of M introduced in (6.4) (with E = E^), we can write 



(8.7) E[e.A«te+V#?>)] =EK-AW(C+V^ 



-E 



2 e-AW(v^-v^) 



We analyse the two terms on the right-hand side separately, beginning with the 
first. Note that (evaluated at the origin) is equal to A° unless one of the edges 
attached to the origin belongs to E^ p \ Moreover, the event that two or more of 
these edges belong to E^ has probability 0(p 2 ). As a consequence, and since 
E[|V0°| 2 ] is bounded uniformly over fj,, we have 



.8) E[£-AW(£ + V0°)] 



EK-A°(£ + V<0)+0(p 2 ). 



e:e=0 



We now turn to the second term in the right-hand side of (8.7), which we decompose 
into 



.9) E 



eG-E< p ',e=0 



■ E 



By Lemma 6.2 and (3.10), it is clear that E[|V0^| 2 ] is bounded uniformly over \x. 

Hence, reasoning as above, we get 

(8.10) 



E 



eeE (p> ,e=0 



= P^E[^.A e (V^-V0°)]+O(p 2 ). 



e:e=0 



We now consider the second term in (8.9). For an edge e such that e ^ 0, the event 
e 6 is independent of the value of A^ (at the origin). Hence, we can rewrite 
the second term in (8.9) as 

P J2 E[£.AW(V^-V0°)], 

e:e^0 
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which itself can be rewritten as 



e:e=0 



where tp^ was introduced in Proposition 7.4. Arguing as above, we see that, on the 
one hand, 

]T EE • A«(V^ - V0°)] = ^ EE ■ A°(V^ - V0°)] + 0(p), 

e:e— e:e— 

while on the other hand, since Dip^ remains bounded in L 2 (Q) by Proposition 7.4, 

EE • A^£>W] = EE • A°Ihl>A + 0{p). 
We have thus shown that 



E 



= pEE • A°Xty M ] " P E E K • ^°(V^ - V0°)] + 0(p 2 ). 

e:e=0 



Combining this with (8.4), (8.5), (8.7), (8.8), (8.9) and (8.10), we thus obtain 
^oU = EE-A°(£ + V^) 



(8.ii) C • 4oU = EE • a°{z + v<o + P <( M ) + o(p), 



where we introduced 

(8.12) al(fx) = £ (EE ' ^ e (C + V^)] - EE ■ A°(£ + V0°)]) + 

e:e=0 

ee • a°d^] - E E ^ ■ A °( y ^ - v^)]. 

e:e=0 

Using the definition of ip^, one can observe that this definition coincides with 
that given in (8.2). Using Theorem 5.2 again, and recalling (8.3), we see that 
EE" A°(£ + V0°)] = £. A£ om £ + o(p). In order to conclude the proof, it thus suffices 
to show that a°(/i) converges to a constant as [i tends to 0. But this is a consequence 
of Propositions 7.2 and 7.4. □ 

9. Approximation by periodization 

The aim of this section is to give alternative characterizations of a\ , based on 
computing correctors on periodizations of the medium. We recall that we write 
B n = {— n, . . . , n} d , and we let B„ = {e € B : e G B n }. We define the periodized 
environment [w]„ G H by letting, for every e € B, ([w]„) e = w e+I , where a; is the 
unique element of (2n + l)Z d such that e + x G B„. By the Lax-Milgram lemma 
and the Poincare inequality (or by basic linear-algebra considerations), for every 
to G O, there exists a unique _B„-periodic function </>° (•, [w] n ) with zero average over 
B n such that for every x E Z d , 

-V* ■A°(Z + V<t>° n )(x > [u,] n ) = 0. 

For the same reason, for every e G B„, there exists a unique -B„-periodic function 
</>*(•, [w] n ) with zero average over B n such that for every x G Z d , 

-v*-A e K + vc)(i,Nn) = o. 

We write E„[/] as a shorthand for E[/([w] n )], whenever this is well-defined. As 
usual, if the function is of the form f(x, [w]„), x G Z d , we interpret E„[/] to mean 

E[/(o,y„)]. 
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Proposition 9.1 (approximation of a\ by periodization). The coefficient a\ ap- 
pearing in Theorem 8. 1 is the limit of a\ (n) as n tends to infinity, where 

(9.1) aUn) = £ ( E » K ' A e (£ + V<0 - E„K • A°(£ + V<0) . 

eSB„ 

Proof. Denning 

«n)= ^(C"C)(0,yn), 

eel,, 

we note that, in close resemblance with (8.12), 

(9.2) al(n) = ]T (E„[£ • A e (£ + V«0 - E„[£ • A°(£ + V<0) + 

e:e=0 

e:e=0 

It thus suffices to prove the following three statements. 

(1) The function uj_ V0° (0, [w] n ) converges in L 2 (fl) to V0°(O, •) as n tends 
to infinity. 

(2) The function oj i-> V<?!>^(0, [w] n ) converges in L 2 (Q) to V</> e (0, •) as n tends 
to infinity. 

(3) The function u i— > D'0„([a;] rl ) converges in £ 2 (H), as n tends to infinity, to 
the function x defined in Proposition 7.4. 

Part (1) is classical (see for instance [CI03, BP04]), but it will be useful to recall 
a proof. We begin by observing that the function 0° is stationary, in the sense 
that for every w e ft and every x G Z d , we have 4> n (x, [w]„) = <^(0, ^ [w] n ). As a 
consequence, we may identify 0° with the function </>° (0, •). We then observe that 

(9.3) E n [D</>° n ■ A°Dr n ] = -E n [Dft ■ , 
and thus, by the Cauchy-Schwarz inequality, 

c_E ti [|,DCI 2 ] < E„ [D4>1 ■ A°Dcj>° n ] < c + |e| E n [\D^ n \ 2 ] 1/2 . 

It follows from this inequality that E„ | 2 ] is bounded uniformly over n, and 

thus that lu H> (0, [w]„) converges weakly in L 2 (H) along a subsequence, to 
some x € i 2 ([l). For notational simplicity, we keep implicit the fact that we now 
consider (Dfi^) only along this subsequence. In order to prove part (1), it suffices 
to show that the weak convergence is actually strong convergence in L 2 (f2), and 
that x = V0°. 

Note first that (9.3) implies that (along the subsequence) 

(9.4) lim E n [Dct>° n ■ A°D<j>° n ] = -E [% • A°$ . 

n— f oo 

Now, take any bounded function / : Q — > E that depends on the value of the 
environment only at a finite number of edges. We have 

E„ [Df ■ A°D<p° n ] = -E„ [Df ■ A $ . 

Passing to the limit (along the subsequence), we get that 

(9.5) E [Df ■ A° x ] = -E [Df ■ A°$ . 

The identity above can then be extended to arbitrary / <E L 2 (Vf). Replacing / by 
</>° in the above identity, and letting n tend to infinity, wc thus learn that 

(9.6) E [ x • A° x ] = -E [x ■ A°d} . 
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To see strong convegence of D<j) n to \ m L 2 (Q), it now suffices to write 

E n [{Dcj> n -x)-A a {Dcj> n -x)} 

= E„ [D<t> n ■ A°D<fi n ] + E„ [x • A° x ] ~ 2E„ [D<t> n ■ A° x ] , 

and observe that the right-hand side tends to zero by (9.4) and (9.6). This ensures 
that x £ Ly, and by (9.5), x satisfies 

-D* ■ A° x = D* ■ A°£. 

By Theorem 7.f , there is a unique such x, and it is V</>° (see Proposition 7.2). 

Part (3) can be obtained in a similar way, noting that the function ip n is stationary, 
in the sense that tp„(x, [oj] n ) — i/j n (0,9 x [w]„). 

For part (2), we cannot argue in the very same way, since A e is not stationary. 
We note instead that, letting cf) n = (p^ — 0° , we have 

at every point (x, [w] n ), with x E Z d and well. We can then simply quote [ALf 1 , 
Lemma A. 2], or argue as in part (1), with the difference that the arguments need to 
be carried out over the physical space 7L d instead of the space of environments. We 
get that (x,ui) (-> lB n V(j> n (x, [w]„) converges in L 2 {Z d x 0) to V^C € L 2 (Z d x 0) 
as n tends to infinity, where satisfies 

(9.7) - V* • A e V^ = V* -C^ + V^ ) (inZ^xO) 

(and the space L 2 (Z d x ft) is defined with respect to the measure obtained as the 
product of the counting measure over Z d and P). This identifies as 

^l(x) = - (V 1 V 2 G e ( a; ,e)) C e (e)(£ + V^°(e)). 

Comparing this with the formula for Wcf> e given in (7.10), we see that we are 
done. □ 

The question arises as to whether the sum involving all (4>n)eem n in the formula 
defining a\{n) can be replaced by a spatial average involving only (<^j) e :e=0' This 
question is important since in practice, one would like to compute as few </>^'s as 
possible. We define 

(9.8) a n (M»)=E ^■A e ^ + ^r n )-^-A°{i + ^<t>l)){^Un). 

x£B n e:e=0 

Proposition 9.2 (approximation of a° by periodization and spatial average). The 
random variable uj i— > a n ([oj] n ) converges in L (O) to the random variable 

^(£ + V0°)-C e (£ + V0 e ) (0,-). 

e:e=0 

Moreover, the coefficient a\ given by Theorem 8. 1 satisfies 

al = V E [(£ + V0°) • C e (£ + V</> £ )] = V lim E„ [(£ + V0°) • C e (£ + V^)] . 

e:e=0 e:e— 

Proof. Let e be such that e = 0. By the definition of 0^ and ^-periodicity, 

We can decompose A e into A° + C e , and then observe that the definition of 0° 
implies that 

xeB n x£B n 
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Since A° is a symmetric matrix, we have obtained that 

(9.9) a n (Mn) = E + V &) ' + V ^)(°' Nn). 

e:e=0 

The convergence announced in the proposition then follows from the fact that 
lu i-> V0° (0, [w]„) and u i-> V<(0, [w]„) converge in L 2 (£i) to V0°(O, •) and V</> e (0, •) 
respectively. The second part of the proposition follows noting that E„[a„] = aj(n), 
see (9.1). □ 



10. First-order expansion around every point 

The aim of this section is to generalize Theorem 8.1 by giving a first-order 
expansion of Ap^t^ as p tends to 0, for every p £ [0, 1]. 

In order to state the result, it is convenient to introduce some new notation. For 
every p £ [0,1] and every e G B, we let A^< e+ ) = A Emu ^ ( reC all that A E was 
defined in (6.1)), <f>^' e = 4>^ u ^ (defined in (6.2)), and so on. Similarly, we 
define A&O = A E ™\^, # ° = 4>f^ {e} , and so on. 

Theorem 10.1 (first-order expansion around every point). Let p £ [0, 1]. There 
exists a\ £ K such that, as p tends to 0, 

(Id) { • A^k = € • AgLt + P af + o( P ). 

Moreover, the coefficient can be defined as the limit of af (p) as fi tends to 0, 
where a®(/i) is given by 

(10.2) af V) - E ' AiV,e+) tt + V ^f' e+) )] " E K • ^ e ~>(£ + V^ e ~>)] 

eGB 

Alternative characterizations of are given in Remark 10.2 below. 

One may attempt to prove this result by defining a coupling of the Bernoulli 
random variables such that if Pi < P2, then C E( p2 \ and then try to adapt 

the proof of Theorem 8.1. This approach leads to a serious difficulty when trying to 
estimate the left-hand side of (6.15). Indeed, this estimate crucially relies on the fact 
that the non-perturbed environment is independent of the events e £ E^f 1 . A naive 
adaptation of this proof to cover the expansion of A^£' as p tends to would thus 

require that the environment be independent of the events e £ E^ +pS> \ 
and of course, this independence does not hold. 

This problem can be overcome with a control of higher moments of the gradients 
of the Green functions than that given by Proposition 3.7. Such results have been 
obtained very recently in the remarkable work [M012], under certain conditions on 
the distribution. We will describe and use these results in the next section, but for 
the problem at hand, a more direct proof can be obtained. 

Proof of Theorem 10.1. It comes out of the above discussion that the decision as to 
which edges are to be perturbed should be taken independently of the environment. 
We thus introduce new independent Bernoulli random variables (6^) e eB of parame- 
ter p G [0, 1], independent of everything else, defined on the same probability space 
(possibly enlarging it). For p < 1 and p £ [0, 1 — p], we further define 

V=Y^ and ^f' v) = (1-^) +6® 
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Recalling that we write for the law of Lie , we see that the law of We p is 



(i - + pv^ = (i - p) [(i - py o) + 



■pi/ 1 ) 



= (1 - p - p)z/ (0) + (p + p)v {1) = v {p+v \ 

In other words, the law of £je is that of u>^P +P \ and as a consequence, the law of 
q (.p,p) ._ (a}f'f)) e£B i s that o£u&*\ 

We can apply Theorem 8.1 for the environment £j(P' p \ seen as a Bernoulli 
perturbation of w®, with the perturbed environment taken from w' 1 ' and the 
perturbation parameter being p. Since the homogenized matrix of that of 



to 



(p+p) 



, we get that as p ^ tends to 0, 



e-4^ = c-41^+pa®+ (p-), 

where is the limit as p tends to of 5^(p) given by 

e£B 

We now see that the quantity 

£ • ^ e+ >(£ + V^ e+ >) - £ • + V0f ) 

is equal to if e £ 25®, and otherwise is equal to 

£ . ^ e+ >(£ + V^f e+ )) - £ • A^Xt + V^ e ")). 
Since this last quantity is independent of the event e € we obtain that 

(10.3) E k- A^ e+ Xt + V<^f' e+) )-£- A®(f + 

= P [e g £®1 E • A^ e+ )(^ + V<?(>jf e+) ) - £ • .4 (p '°(£ + Wjf'°) 

since P [e £ £■<»] = 1 — p, we obtain the desired result for p ^ (i.e. for the right 

derivative of £ • A^ m £). 

A similar reasoning can be followed if p € (0, 1] and p € [— p, 0], by letting 

p==? and ,;»> = (l-^<» + ^<°> 

The reasoning leads to the same formula for the left derivative, so the proof is 
complete. □ 

Remark 10.2. From the construction in the proof of Theorem 10.1 and Proposi- 
tion 9.1, one can give an alternative characterization of the coefficient a® appearing 
in Theorem 10.1. Indeed, it is the limit of a\(n) as n tends to infinity, for 

a x {n) = £ (En [e-^' e+) (e + V^ e+ ')] -E„ [f • + V# e ">)]) , 

eGB„ 

where (0 ,e ' and </>l P ' e ' are defined as 0°, but with respect to A^ ,e+ ) and A^ ,e ) 
respectively instead of A° . 

Similarly, using Proposition 9.2, we can obtain yet another characterization of 
a i , namely 



(10.4) 
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where V<^>' e )(0,-) and V<^' e+ )(0, •) are the L 2 (f])-limits of V^' e } (0,-) and 
^4>n' e '(0, ■) respectively, as n tends to infinity. Moreover, 

(10.5) = ^ lim E„ [(£ + V^< e ~)) • G e (£ + V^' e+) ) 

* — » n— foo L 

e:e=0 

11. Refined control of the error term 

The aim of this section is to obtain a sharper control of the error term o{p) 
appearing in (10.1) of Theorem 10.1. This generalization will be achieved only at 
the price of an additional assumption. 

As is clear from the proof of Theorem 10.1, in order to achieve this, it is in fact 
sufficient to control the error term in Theorem 8.1. The fact that Theorem 8.1 gives 
only a weak control of the error term lies in Theorem 6.1. The weak point in the 

argument leading to Theorem 6.1 is in the estimation of the left-hand side of (6.15), 

(p) 

where we crudely estimated the second mixed derivatives of G/f using the pointwise 
inequality (3.10) on the gradient of Gp f . Although Proposition 3.7 shows that the 
L 1 norm of the second mixed derivatives of Gjj have same spatial decay as the 
homogeneous case, this is of no use when estimating the left-hand side of (6.15) 
because we have no control on the dependence between the second derivatives of 
G^ p) and the events e € . 

This problem can however be overcome with a control of higher moments of the 
second derivatives of the Green function. Such results have been obtained very 
recently in the remarkable work [M012], under certain conditions that we now 
describe, starting with some definitions. 

For every bounded measurable function / : [c_, c + ] — > R + , we define the entropy 
of / with respect to the law of u> e as 

Ent(/) =E[/(w e )log/(w e )] -E[/(w e )]logE[/(w e )]. 

We say that the law of cj e satisfies a continuous log-Sobolev inequality with 
constant p > if for every continuously differentiable function / : [c_, c + ] — > K, we 
have 

Ent(/ 2 ) < -E[(/'(c e )) 2 ]. 
P 

We say that the law of uj e satisfies a discrete log-Sobolev inequality with constant p > 
if for every bounded measurable function / : [c_,c + ] — > K, we have 

Ent(/ 2 K -E[(/(c e )-E[/K)]) 2 ]. 
P 

As is well-known, a very convenient feature of the log-Sobolev inequalities is that 
they can be tensorized to yield inequalities for functions defined on the product 
space ft (see [GZ03, Theorem 4.4]). Such inequalities form one of the building blocks 
of the approach of [MO 12]. We only quote the result of [MO 12] that concerns the 
second mixed derivative of the Green function, but similar results were obtained for 
the first derivative as well. 

Theorem 11.1 ([M012]). Assume that the law of ui e satisfies either the continuous 
or the discrete log-Sobolev inequality with constant p > 0. For every q ^ 1, there 
exists a finite C = G(d, c_, c+, p, q) such that for every leZ" 1 and p, ^ 0, 

G 



E[\V 1 V 2 G ll (0,x)\i] 1/q 



(i v \x\y 

Remark 11.2. Striclty speaking, [M012, Theorem 1] corresponds to this result 
for p = 0. The proof can however be adapted with minor modifications to yield 
Theorem 11.1. Let us describe briefly how. [M012, Lemma 1] needs no change. 
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[M012, Lemma 3] consists of two parts. The proof of the first part can be easily 
adapted to yield the result with G replaced by G M (actually, it provides an alternative 
route to the periodization argument used there), while the second part was proved 
for G M in the first place (that is, in [GOll]). The first part of [M012, Lemma 2] 
is all what is needed for our purpose, and it remains true with G replaced by G M 
(uniformly over /i). Indeed, the formulas for the derivatives of the Green function 
with respect to ui e appearing in steps 1 to 4 of the proof, as e.g. [M012, (39)-(41)], 
remain valid provided G is replaced by G^ everywhere. Step 5 is a consequence 
of [M012, Lemma 3], while step 6 follows from previous steps, and concludes the 
proof of the part of the lemma we are interested in. Finally, the proof of [M012, 
Theorem 1] follows from the results we just reviewed (again providing an alternative 
route to the periodization argument proposed there). 

The results we will prove in this section are conditional on an assumption that 
we now introduce, which is related to the conclusion of Theorem 11.1. For every 
p € [0, 1] and every q ^ 1, let G (0, +oo] be defined by 

(11.1) C<>)= sup (lV|z|) d E[|V 1 V 2 G(f ) (0,2;)H 1/ \ 

Let p £ [0, 1]. We say that assumption (H® ) holds if for every q ^ 1, there exists 
a neighbourhood of p in [0, 1] such that is bounded uniformly over p in this 
neighbourhood. 

Theorem 11.3 (sharp control of the error). Let p <E [0, 1] and <G K be given by 
Theorem 10.1. If assumption (H® ) holds, then for every r\ > 0, 

(ii-2) e • = t ■ 4*U + p «?° + °(^"). 

Before we proceed to the proof of Theorem 11.3, we comment on the relations 
between assumption (H®) and Theorem 11.1. We see from Theorem 11.1 that in 
order to ensure that (H®) holds, it suffices to show that the law of u?' satisfies a 
continuous or discrete log-Sobolev inequality with constant p( p \ and that is 
uniformly bounded away from in a neighbourhood of p. 

For many distributions of interest, as for instance all those with a finite support, 
the continuous log-Sobolev inequality fails to hold. We thus turn our attention to 
the discrete version of the inequality. The set of distributions satisfying a discrete 
log-Sobolev inequality is well-understood. Indeed, the law of u> e satisfies a discrete 
log-Sobolev inequality if and only if it has finite support. In this case, the optimal 
constant in the log-Sobolev inequality is explicit, and given by 

1 - 27T* 

( ' ]og[(l-7rO/7r*]' 

where 7T* = minP[cj e = x], and the minimum is taken over all x ranging in the 
(finite) support of the law of u e [DS96, Theorem A.l] (here, we assumed that the 
law of u> e is not concentrated on a single point, in which case 7T* ^ 1/2 ; the quantity 
in (11.3) can be extended by continuity to 71% = 1/2 and takes the value 1/2 there). 

The constant in (11.3) is bounded away from as soon as 71% itself remains 
bounded away from 0. By Theorem 11.1, this shows that assumption (H®) holds 
as soon as the laws of uj e and wi 1 ' both have finite support and p {0, 1}. If the 
laws of uj e and ui 1 ' both have finite support, then assumption (H°) holds as soon as 
the support of lo^ is included in that of w e (we write (H°) for (H^)). The same 
statement with w e and uj^ 1 interchanged holds concerning assumption (H^). 

Note for instance that the argument above cannot show that assumption (H°) 
holds in the case where the environments are (distinct) constants! Yet, the idea 
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that assumption (H°) really breaks down in this case is highly doubtful. It is thus 
very likely that the problem we are encountering here is an artefact of the proof 
techniques. In fact, I believe that assumption holds for every p without any 

further assumptions than those introduced in section 2. 

We now turn to the proof of Theorem 11.3. As was said before, it suffices to 
prove it for p = 0, and the key step for doing so is to prove the following result. 

Theorem 11.4 (sharp linear approximation of the corrector). If assumption (H°) 
holds, then for every r\ > 0, there exist constants c, 7 > such that for every 
fx e (0,1/2] andpe [0,1], 



E 



We begin the proof with a result whose purpose is to replace Lemma 6.5. 

Lemma 11.5 (sharp averaged control of the perturbed Green function). There 
exists a constant c < +00 ( depending only on the ellipticity constants and on the 
dimension) such that for every q ^ 1, fi G (0, 1/2] and e,e' G B, 

E[|V 1 V 2 G^(e',e)r] 1/9 S =: C Q ' 



(IV |e-e'|)° 



Proof. Recall from the proof of Lemma 6.5 that — — G°. Since by definition 



of C°, 



1/9 / S 



(11.4) E[\V 1 V 2 G;(e',e)\^ L/q 
it suffices to prove that 



(IV |e-e'|) d ' 



E 



|ViV 2 G'(e',e)|« 



1/9 



^ c 



(1 V \e-e'\) d ' 



Recall from (6.10) that 

ViVaS^e) = G e (e) (V 1 V 2 G*(e ) e)) (V 1 V 3 GP(e' J e)) 
Since, by (3.10), 

|G e (e)V 1 V 2 G; t (e,e)| 
is bounded uniformly over /x and w, the result follows from (11.4). 



□ 



Proof of Theorem 11. 4. We assume that (H°) holds. Let q > 3. By the assumption, 
the numbers defined in (11.1) are bounded by some constant C q uniformly over 
p ^ in a neighbourhood of 0. From now on, we assume that p is sufficiently small 
to lie in this neighbourhood. 

For the beginning of the proof, we can follow the same reasoning as in the proof 
of Theorem 6.1, up to the point when we arrive at the estimation of the left-hand 
side of (6.15), that is, 



E 



£ |ViV 2 G(f)(0,e')| |V 1 V 2 G«(e',e)| |£ + V<£°(e)| 



where we recall that E^ = E^ n B^--,, and 7 > 1/2. We rewrite slightly this 
expectation as 

(11.5) £ E 



IVrVsG^e')! iVxVaG^e)! |£ + V0°(e) 
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Let £ = 1 — 3/q. We apply Holder's inequality with exponents (£ , q, q, q) to bound 
each summand in (11.5) by 

(11.6) P[e,e'e#] ( E 

E [ |V 1 V 2 G*(e',e)|*] 1/9 E [|£ + V0°(e)|*] 1/9 . 

Since we only consider summands for which e^e',we have 

P[e,e' € £ (p) ] c =p 2C . 

By our assumption and Lemma 11.5, the product of the second and third terms in 
(11.6) is bounded by 

(IV \e/\) d (IV |e-e'|) d ' 

while, by Theorem 5.2, the last term in (11.6) is bounded by some power of log(/x _1 ), 
say log r (/i _1 ). We thus obtain that the sum in (11.5) is bounded by 

By comparing this sum to an integral, we get that the whole expression above is 
bounded, up to a constant, by p 2 ^/i~ 7 '', for any -q > 0. Recalling that ( = 1 — 3/q 
and that q can be chosen arbitrarily large, this concludes the proof. □ 

Proof of Theorem 11.3. As noted before, it suffices to prove the theorem for p = 0. 
The proof follows closely that of Theorem 8.1. Let r\ > 0. Instead of (8.3), we fix 
fi = p 2 . With this choice, Theorem 5.2 implies that the left-hand side of (8.4) is 
(Dip 2 -* 1 ), while by Theorem 11.4, the left-hand side of (8.5) is 0(p 2 - 3r >). We can 
then follow the proof of Theorem 8.1 without any change until (8.11), whose error 
term is now 0(p 2 ~ 3ri ) instead of o(p). The proof is then concluded in the same way, 
noting that Theorem 5.2 ensures that E[£ ■ A°(£ + V<£°)] = £ ■ A° lom ^ + 0{p 2 ^). □ 

Acknowledgements. I would like to warmly thank Vincent Beffara, Noam Berger, 
Antoine Gloria, Claude Le Bris and Frederic Legoll for stimulating discussions about 
this problem. 

References 

[A113] Y. Almog. Averaging of dilute random media: a rigorous proof of the Clausius— Mossotti 

formula. Arch. Ration. Mech. Anal. 207 (3), 785-812 (2013). 
[ALlOa] A. Anantharaman, C. Le Bris. Homogeneisation d'un materiau pcriodique faiblement 

perturbe aleatoirement. C. R. Math. Acad. Set. Paris 348 (9-10), 529-534 (2010). 
[ALlOb] A. Anantharaman, C. Le Bris. Elements of mathematical foundations for a numerical 

approach for weakly random homogenization problems. Preprint, arXiv: 1005.3922 (2010). 
[ALII] A. Anantharaman, C. Le Bris. A numerical approach related to defect-type theories for 

some weakly random problems in homogenization. Multiscale Model. Sirnul. 9 (2), 513—544 

(2011). 

[Be05] J. Beltran. Regularity of diffusion coefficient for nearest neighbor asymmetric simple 
exclusion on Z. Stochastic Process. Appl. 115 (9), 1451-1474 (2005). 

[BM01] L. Berlyand, V. Mityushev. Generalized Clausius-Mossotti formula for random composite 
with circular fibers. J. Statist. Phys. 102 (1-2), 115-145 (2001). 

[Be02] C. Bernardin. Regularity of the diffusion coefficient for lattice gas reversible under Bernoulli 
measures. Stochastic Process. Appl. 101 (1), 43-68 (2002). 

[BLL07] X. Blanc, C. Le Bris, P.-L. Lions. Stochastic homogenization and random lattices. J. 
Math. Pares Appl. 88 34-63 (2007). 

[BP04] A. Bourgeat, A. Piatnitski. Approximations of effective coefficients in stochastic homoge- 
nization. Ann. Inst. H. Poincare Probab. Statist. 40 (2), 153-165 (2004). 



|ViV a G<?>(0,e')|« 



EXPANSION OF HOMOGENIZED COEFFICIENTS 



31 



[CI03] P. Caputo, D. Ioffe. Finite volume approximation of the effective diffusion matrix: the case 

of independent bond disorder. Ann. Inst. H. Poincare Probab. Statist. 39 (3), 505-525 (2003). 
[CNOOa] J. Conlon, A. Naddaf. Green's functions for elliptic and parabolic equations with random 

coefficients. New York J. Math. 6, 153-225 (2000). 
[CK08] M. Cudna. T. Komorowski. A finite dimensional approximation of the effective diffusivity 

for a symmetric random walk in a random environment. J. Comput. Appl. Math. 213 (1), 

186-204 (2008). 

[De99] T. Delmotte. Parabolic Harnack inequality and estimates of Markov chains on graphs. Rev. 

Mat. Iberoamericana 15 (1), 181-232 (1999). 
[DD05] T. Delmotte, J.-D. Deuschel. On estimating the derivatives of symmetric diffusions in 

stationary random environment, with applications to V0 interface model. Probab. Theory 

Related Fields 133 (3), 358-390 (2005). 
[DS96] P. Diaconis, L. Saloff-Coste. Logarithmic Sobolev inequalities for finite Markov chains. 

Ann. Appl. Probab. 6 (3), 695-750 (1996). 
[EGMN12] A.-C. Egloffe, A. Gloria, J.-C. Mourrat, T.N. Nguyen. Random walk in random 

environment, corrector equation, and homogenized coefficients: from theory to numerics, back 

and forth. Preprint, arXiv:1211.1834. 
[GOll] A. Gloria, F. Otto. An optimal variance estimate in stochastic homogenization of discrete 

elliptic equations. Ann. Probab. 39 (3), 779-856 (2011). 
[G012] A. Gloria, F. Otto. An optimal error estimate in stochastic homogenization of discrete 

elliptic equations. Ann. Appl. Probab. 22 (1), 1-28 (2012). 
[GZ03] A. Guionnet, B. Zegarlinski. Lectures on logarithmic Sobolev inequalities. Seminaire de 

probabilities XXXVI 1-134, Lecture notes in math. 1801, Springer (2003). 
[Hu09] B. Hughes. Conduction and diffusion in percolating systems. Encyclopedia of complexity 

and systems science 1395-1424 (2009). 
[KV86] C. Kipnis, S.R.S. Varadhan. Central limit theorem for additive functionals of reversible 

Markov processes and applications to simple exclusions, Comm. Math. Phys. 104, 1-19 (1986). 
[Ko89] S.M. Kozlov. Geometric aspects of averaging. (Russian) Uspekhi Mat. Nauk 44 (2), 79-120 

(1989). Transl. in Russian Math. Surveys 44 (2), 91-144 (1989). 
[KO05] T. Komorowski, S. Olla. Einstein relation for random walks in random environments. 

Stochastic Process. Appl. 115 (8), 1279-1301 (2005). 
[Kii83] R. Kiinnemann. The diffusion limit for reversible jump processes on Z d with ergodic random 

bond conductivities. Comm. Math. Phys. 90 (1), 27-68 (1983). 
[LOV01] C. Landim, S. Olla, S.R.S. Varadhan. Symmetric simple exclusion process: regularity of 

the self-diffusion coefficient. Comm. Math. Phys. 224 (1), 307-321 (2001). 
[LOV04] C. Landim, S. Olla, S.R.S. Varadhan. On viscosity and fluctuation-dissipation in exclusion 

processes. J. Statist. Phys. 115 (1-2), 323-363 (2004). 
[M012] D. Marahrens, F. Otto. Annealed estimates on the Green's function. Preprint, available 

at http://www.mis.mpg.de/de/publications/preprints/2012/2012-69.htinl (2012). 
[Ma] J.C. Maxwell. Medium in which small spheres are uniformly disseminated. A treatise on 

electricity and magnetism, 3d ed., part II, chapter IX, article 314. Clarendon press (1891). 
[Moll] J.-C. Mourrat. Variance decay for functionals of the environment viewed by the particle. 

Ann. Inst. H. Poincare Probab. Statist. 47 (1), 294-327 (2011). 
[Na05] Y. Nagahata. Regularity of the diffusion coefficient matrix for the lattice gas with energy. 

Ann. Inst. H. Poincare Probab. Statist. 41 (1), 45-67 (2005). 
[Na06] Y. Nagahata. Regularity of the diffusion coefficient matrix for generalized exclusion process. 

Stochastic Process. Appl. 116 (6), 957-982 (2006). 
[Na07] Y. Nagahata. Regularity of the diffusion coefficient matrix for lattice gas reversible under 

Gibbs measures with mixing condition. Comm. Math. Phys. 273 (3), 637-650 (2007). 
[Pa95] G.C. Papanicolaou. Diffusion in random media. Surveys in applied mathematics 205—255 

(J.B. Keller et al. Eds), Plenum Press (1995). 
[PV81] G.C. Papanicolaou, S.R.S. Varadhan. Boundary value problems with rapidly oscillating 

random coefficients. Random fields (Esztergom, 1979) 835-873, Colloq. Math. Soc. Janos 

Bolyai 27, North-Holland (1981). 
[SC] L. Saloff-Coste. Lectures on finite Markov chains. Lectures on probability theory and statistics 

(Saint-Flour 1996) 301-413, Springer (1997). 
[SW01] S. Smirnov, W. Werner. Critical exponents for two-dimensional percolation. Math. Res. 

Lett. 8 (5-6), 729-744 (2001). 
[SZ97] D.W. Stroock, W. Zheng. Markov chain approximations to symmetric diffusions. Ann. Inst. 

H. Poincare Probab. Statist. 33 (5), 619-649 (1997). 
[St] J.W. Strutt (3d Baron Rayleigh). On the influence of obstacles arranged in rectangular order 

upon the properties of a medium. Philos. mag. 34, 481-502 (1892). 



32 



JEAN-CHRISTOPHE MOURRAT 



[Su05] M. Sued. Regularity properties of the diffusion coefficient for a mean zero exclusion process. 

Ann. Inst. H. Poincare Probab. Statist. 41 (1), 1-33 (2005). 
[ZB77] M. Zuzovsky, H. Brenner. Effective conductivities of composite materials composed of 

cubic arrangements of spherical particles embedded in an isotropic matrix. Journal of applied 

mathematics and physics (ZAMP) 28, 979-992 (1977). 

EPFL, INSTITUT DE M ATHEM ATIQUES , STATION 8, 1015 LAUSANNE, SWITZERLAND 



